Influence of trigonal warping on interference effects in bilayer graphene 
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Bilayer graphene (two coupled graphitic monolayers arranged according to Bernal stacking) is a 
two-dimensional gapless semiconductor with a peculiar electronic spectrum different from the Dirac 
spectrum in the monolayer material. In particular, the electronic Fermi line in each of its valleys 
has a strong p — > — p asymmetry due to trigonal warping, which suppresses the weak localization 
effect. We show that weak localization in bilayer graphene may be present only in devices with 
pronounced intervalley scattering, and we evaluate the corresponding magnetoresistance. 

PACS numbers: 73.23.-b, 72.15.Rn, 73.43. Qt, 81.05.Uw 
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Weak localization (WL), which universally occurs in a 
broad range of disordered conductors [l|, Q is caused by 
the constructive interference of electron waves circling 
the same closed path in opposite directions. It is sen- 
sitive to the presence of an external magnetic field as 
manifested in a negative magnetoresistance (MR) at low 
temperatures. 

Usually qualitative features of WL do not depend on 
the detail of the electronic band structure and crystalline 
symmetry, with the exception of spin-orbit coupling 0, 
3] . In gapful multi- valley semiconductors only the size of 
the WL effect may depend on the number of valleys and 
the strength of intervalley scattering 0, d, [f| . The low- 
field MR, Ap(B) = p(B) — p(0), in a two dimensional 
electron gas or a thin metallic film [3, 0, 0, m the 
absence of spin-orbit coupling is characterized by 
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Here F{z) = lnz + ^(5 + j), 4>{z) is the digamma func- 
tion, t v is the coherence time, D is the diffusion coeffi- 
cient, and the integer factor sg depends on whether or not 
states in n v valleys are mixed by disorder. This factor is 
controlled by the ratio 6 = ti/t^ between the intervalley 
scattering time t\ and the coherence time t v . In mate- 
rials such as Mg, ZnO, Si, Ge (listed in Table I) where 
each of the Fermi surface pockets is p — > p symmet- 
ric, intervalley scattering reduces the size of the WL MR 
from that described by = 2n v when 9 — Ti/r„ 3> 1 to 
s = 2 for 9 < 1. 

A more interesting scenario develops in a multi-valley 
semimetal, where the localization properties can be in- 
fluenced by the absence of p — > p symmetry of the 
electronic dispersion within a single valley. Graphene 
.10, [Til HH in a graphene-based transistor [13, [lj, [l5j 
represents an example of such a system. In this Letter, 
we demonstrate how the asymmetry in the shape of the 
Fermi surface in each of its two valleys determines the 
observable WL behavior. In contrast with conventional 
materials (see Table I), WL MR in bilayer graphene is 
increased by the intervalley scattering, from Soo = at 
9 — > 00 to s = 2 at 9 < 1. 

Bilayer graphene consists of two coupled graphitic 



monolayers arranged according to Bernal stacking, see 
[l2| for details of the lattice configuration. Its unit cell 
contains sites A, B and A, B. Sites A, B and A, B belong 
to the honeycomb lattice in the bottom and top layers, 
respectively, with sites B being exactly below A. The 
hexagonal Brillouin zone of the bilayer has two inequiv- 
alent corners K + and K_ The four branches of its 
electronic spectrum [17] form one pair split by about ±71 
(the interlayer coupling) and two low-energy branches 
(formed by states based upon sublattices A and B) which 
are degenerate at K + and K_. The low-energy branch 
can be described [T3] using a Hamiltonian, 
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which acts in the space of four-component wave func- 
tions $ = [</>k + ,a,0 k+ ,b>0k_,s>0k_,ji]- Here, is 

an electron amplitude on the sublattice a = A, B in the 
valley £ = K + ,K_, tJ x y-z and H x , y .z a- re Pauli matrices 
acting in sublattice and valley space, respectively p^ |. 
In slightly doped graphene, two disconnected Fermi lines 
surround the corners of the Brillouin zone. 

The first term in Eq. |T]) is the leading contribution in 
the nearest nei ghb or approximation of the bilayer tight 
binding model [12]. It includes intralayer hopping A ^ 
B and A ^5 B (that leads to the dispersion e = ±pv 
near K± in a monolayer) and interlayer A ^ B hopping, 
and it yields a parabolic spectrum e = ±p 2 /2m with 
m = ji/2v 2 . It has been noticed [l2j that quasiparticles 



TABLE I: WL factor sg in conductors with a multi- valley con- 
duction band and negligible spin-orbit coupling. The factor 
se is specified for two limiting cases, no inter-valley scattering 
6 = Ti/r v — > 00, and for fast inter- valley scattering 6 — > 0. 
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Mg films [8], ZnO wells [9] 
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Si MOSFETs [4, 5] 
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Si/SiGe wells [6j 
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described by this term are chiral [18j : their plane wave 
states are eigenstates of an operator crn 2 with crn 2 = — 1 
for electrons in the conduction band, and aw-i = 1 for 
the valence band, where n 2 (p) = (cos(2<p),sin(2<p)) for 
p = (pcos</?,psin<£>). For an electron in the conduction 
band, the plane wave state is 
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where | T)k+. p = [1,0,0,0], U) K+lP = [0,1,0,0] and 
I T)k_, p = [0,0,1,0], | T>K_.p - [0,0,0,1], and the fac- 
tors e ±4V take into account the chirality. 

The term h w in Eq. ([T]) originates from a weak di- 
rect A ^=5 B interlayer coupling. This one gives rise to 
p — > — p asymmetry in the dispersion of electrons. It also 
leads to a Lifshitz transition in the shape of the Fermi 
line of the 2D electrons which takes place when ep ~ el = 
\li{vz/v) 2 thus n e > n L w v%-yf /(2nh 2 v 4 ) ~ 10 n cm- 2 
(using v 3 /v ~ 0.1, Ref. [12]). As long as n e < n h , the 
Fermi line in each valley is split into four pockets [12[ . For 
n e > rtL, h w can be treated as a perturbation leading to a 
trigonal deformation of a single-connected Fermi line and 
asymmetry of the electron dispersion inside each valley 
illustrated in Fig. 1: e(K+, p) ^ e(K±, — p), though, due 
to time-reversal symmetry [19|, e(K±,p) = e(K T ,— p). 
Existing graphene structures are strongly affected by 
charges trapped in the underlying substrate or on its sur- 
face. Such disorder (as well as the inhomogeneity of the 
electron density) obscures the intricate details of the elec- 
tronic Fermi surface at the lowest energies, e <C £l- How- 
ever the interplay between the two terms in H2L result- 
ing in intravalley asymmetry of the electronic dispersion 
manifests itself in the WL behavior. 

The WL correction to conductivity in disordered con- 



K / 



\ K 



+ 



e(K,p)=e(K + ,-p) 



E(K,p)# 8 (K,-p) 




s(K + ,p)^e(K + ,-p) 



\ / 



FIG. 1: Fermi lines (solid lines) in the vicinity of two inequiv- 
alent valleys K+ and K_ of the hexagonal Brillouin zone 
(dashed line). Trigonal warping produces asymmetry of the 
dispersion at each valley e(K±,p) 7^ e(K±,— p), where mo- 
mentum p is determined with respect to the center of the 
valley, but the effects of warping in the valleys have opposite 
signs, e(K±,p) = e(K T ,— p). The bilayer lattice configura- 
tion is described in detail in 1121. 



ductors is a result of the interference of electrons prop- 
agating around closed loops in opposite directions [ij. 
Such interference is constructive in materials without 
spin-orbit coupling, since electrons acquire exactly the 
same phase when travelling along two time-reversed 
paths. It is usually described [1[ in terms of the particle- 
particle correlation function, Cooperon. To evaluate WL 

in bilayer graphene, we use a Cooperon matrix C^Mofai 
where subscripts label the sublattice state of incoming 
a(3 and outgoing a'/3' pairs of electrons and superscripts 
describe the valley state of incoming and outgoing 
pairs. Following the example of Cooperons for a spin | 
[2| , we classify Cooperons as singlets and triplets in terms 
of sublattice and valley indices Cg 1 ^ 2 . For example, 
M = is a "valley-singlet" , M = x,y, z are three "valley- 
triplet" components; S = is a " sublattice-singlet" and 
S = x,y,z are " sublattice-triplet" components. 

With regards to the sublattice composition of Cooper- 
ons in disordered bilayer graphene, only modes C^' 1 
(with S = z [20]) are relevant. Indeed, a correlator 
C ~ ^K+.p^K-^p describing two plane waves, ^K+.p 
and ^K-.-p, propagating in opposite directions along 
a ballistic segment of a closed trajectory contains only 
sublattice-triplet terms, 
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and the terms corresponding to C% y disappear after av- 
eraging over the direction of momentum, (e ±2lv ) ip = 0. 

With regards to the valley composition of relevant 
Cooperon modes, the symmetry of the electronic dis- 
persion within each valley plays a pivotal role. For 
a conventional metal @ or semiconductor 0, H, @, H|, 
two phases and $q acquired by an electron while 
propagating in a clockwise and anti-clockwise direction 
along the same loop are exactly equal, so that the in- 
terference enhances backscattering and leads to WL [l[. 
Any closed trajectory is a combination of ballistic in- 
tervals characterized by momenta ±pj (for two direc- 
tions) and time of flight tj. The asymmetry of the elec- 
tron dispersion in bilayer graphene, due to h w , gener- 
ates a phase difference S = 1)^ ~ i)q — Y^^j) where 
Sj ~ [e(Pj) — e(— Pj)]*j = 2/i w (pj)£j. Since Sj are ran- 
dom uncorrelated, the mean square of S accumulated 
over the time interval t — ^ . tj can be estimated as 

(S 2 ) ~ ((tjh w (pj)) 2 )t/T, where r is the transport time 
and I = vft is the mean free path [l^]. This determines 
the relaxation rate, 



^T(Tr^(p)) v =nn L l 2 T- 
1 



2ti 
T 



im L l 2 < 1 
■ktilI 2 > 1 



(2) 

which suppresses the intravalley Cooperons Cf, C% . We 
estimate that for the recently studied bilayers [ll| with 
n e = 2.5 x 10 12 cm~ 2 , r w ~ r and I ~ 0.1/zm. A 
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similar situation occurs in bilayer structures studied by 
R. Gorbachev et al. [2l[. Due to time-reversal symme- 
try [e(K±,p) = e(K T ,— p)] signs of warping in the 
valleys K + and K_ are opposite, so that it does not sup- 
press the intervalley Cooperons C z and C z . 

The above qualitative analysis indicates that strong 
trigonal warping should be associated with the suppres- 
sion of WL in bilayer graphene when electrons do not 
change their valley state upon scattering. To give a 
quantitative description of this prediction we used the 
diagrammatic technique and related the WL correction 
to conductivity to the surviving Cooperon modes Q, 
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The corresponding scattering rates 
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where C (r, r) are the Cooperon propagators taken at co- 
inciding coordinates. For completeness, in Eq. ([3]) we 
have retained the intravalley Cooperons C^' v , though 
they are strongly suppressed by trigonal warping. Follow- 
ing their suppression, the WL correction is determined 
by the intervalley modes C z and C z but, in the absence 
of intervalley scattering, the contributions of C z and C" 
are equal in magnitude, so that they cancel. Intervalley 
scattering due to atomically sharp scatterers breaks this 
exact cancellation and partially restores the WL effect. 

A general form of time-reversal-symmetric [l9j dis- 
order in bilayer graphene can be represented follow- 
ing the method used in the monolayer graphene studies 
[22, Hil, HH • That is, we determine irreducible representa- 
tions of the symmetry group of the bilayer crystal formed 
by 4 x 4 matrices acting in the basis of states | T)k + ,p, 
I I)k +:P , I 4)k_, p , I T)k_, p , and, then, use them to con- 
struct the general t — > —t symmetric disorder u(r): 

u = u(r) +y^E s A;it 3 ;(r); s,l = x,y,z; (4) 
Z x = U z <g> a x , Y<y = Il z (g) <jy, Y> z = n <E> a z , 
A x = Il x ® <r z , Ay = n y ® a z , A z = IL Z (g) do- 
Here, E = (E X) E 3/ ,E 2 ) and A = (A x ,A y ,A z ) 
form two mutually commuting algebras [E Sl ,E S2 ] = 
2ie s ^- s ^ S3 , [A h ,A h ] = 2ie l ^A h , and [E, A] = 0. 
Both E s and A; invert sign upon time reversal [l9j |. 
whereas their products remain invariant under t —* —t 
transformation. The first term in u represents the po- 
tential of remote charges which is made short range by 
screening by 2D electrons (u(y)u(y 1 )) ~ u 2 S(r — r'). It 
affects electrons in both layers equally and manifests it- 
self through the scattering rate t^ 1 = tx^u 2 jh 2 where 
7 = m Up jlT ^ 2 ) is the density of states per spin in each 
valley [7|. Disorder u zz generates a random difference 
between energies on the A and B sites (bottom and top 
layers). The two terms containing u xz and u yz origi- 
nate from fluctuations in hopping and they scatter elec- 
trons within each valley, whereas the other terms in u 
take into account intervalley scattering. For simplicity, 
we assume that different types of disorder are uncorre- 
cted, (u s i(r)u s >i'(r')) = u 2 sl 5 ss '6u>6(r - r'), and, on av- 



t s1 x = nju 2 sl /h, where tJ = r^ 1 = t s I and tJ 



can be combined into the intervalley scatter- 
2t~} and the intravalley rate 
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t^ 1 = t zz both of which lead to an additional sup- 
pression of intravalley modes. Intervalley scattering also 
leads to the relaxation of C z although it does not affect 
the valley-symmetric mode C z . All the scattering mech- 
anisms limit the transport time r _1 = Tq 1 + . 

Two low-gap modes C z and C* obey the following 
equations, 
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C°(r,r') = <Kr-r'), 
<7|(r,r')=<5(r-r'), (5) 



where we included an external magnetic field, B = rotA 
in P = (iV + % A), and inelastic decoherence, r~ 1 (T). 

Equations (|3I5|) yield the zero field WL correction to 
the resistivity and the WL MR, 
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erage, isotropic in the 
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where B vi — he/ (ADeT v:i ). Equation © gives a com- 
plete description of the crossover between two extreme 
regimes mentioned at the beginning [7|. It also in- 
cludes small contributions of the suppressed intraval- 
ley Cooperons, So = [2e 2 p/(irh)} ln(r v T*/[r(T* + t v )]) 
and 5(B) = -[2e 2 p/(irh)}F[B/(B v + B*)], where r' 1 = 
t^ 1 + 2t~ 1 +T; - and B* = fic/(4Z3er*). This permits us 
to account for a possible difference between the warping 
time t w and the transport time r. According to Eq. © 
WL MR in bilayer graphene sheet disappears as soon as 
7] exceeds t v , whereas in structures with t v > 71, the 
result Eq. © predicts the WL behaviour, as observed 
in [2l|. Such WL MR is saturated at a magnetic field 
determined by the intervalley scattering time, instead of 
the transport time as in usual conductors, which provides 
the possibility to measure n directly. 

Also, it is interesting to consider a small device of bi- 
layer graphene with both sizes, length L and width L±, 
less than the electron coherence length ^JDt v . Since the 
edge of graphene scatters electrons between K + and K 
valleys, in a wire ("ribbon") with L± < L the sample 
width starts playing the role of the intervalley scatter- 
ing length, even though intervalley scattering in the bulk 
of the material may be irrelevant, L ; = y 1 Dt\ 3> L. In 
this case, in a wire with conductance G = G c iassicai + 
Gwl + 5G both WL and universal conductance fluctua- 
tions (UCFs) [Ifl , would be described by the usual results 
established for the orthogonal symmetry class in disor- 
dered systems [25|, where UCFs and the WL parts are 

of the same order, G WL f£, (5G 2 ) 1 ' 2 - ^jr- 

In contrast, in a short and wide strip with L± > L and 
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no intervalley scattering (either off defects in graphene 
or contacts), WL in most of the cross-section is sup- 
pressed by warping, whereas the UCFs (which are im- 
mune to such symmetry breaking) are not. Therefore, 

in a strip G wh ~ (SG 2 ) 1/2 ~ 4 VtFx> whicn not 

only reflects a larger variance of UCFs in broad samples 
(L± ^> L), but also the fact that in a valley-degenerate 
system two valleys give coherent contributions towards 
the observable conductance fluctuations. 

In conclusion, we have shown that p — > p asymme- 
try of the electron dispersion in each valley of bilayer 
graphene leads to unusual (for conventional disordered 
conductors) behavior of interference effects in electronic 
transport. Without intervalley scattering, trigonal warp- 
ing of the electron dispersion near the center of each val- 
ley leads to a suppression of weak localization while in- 
tervalley scattering restores it. This behavior is expected 
in experimentally available bilayer graphene structures 
[Til [2l| . for which r w ~ t. In contrast, the univer- 
sal conductance fluctuations are not reduced by trigonal 
warping and may be even stronger in a system without 
intervalley scattering than in systems with short r;: a 



completely opposite behavior to the WL effect. Such 
behavior of bila yer grap hene is similar to that of mono- 
layer graphene [22, l2q |. despite the fact that electrons 
in these two materials have different chiralities and can 
be attributed different Berry phases: n in monolayers, 
2tt in bilayers [H, 13 ■ More generally a suppressed 
weak localization magnetoresistance and its sensitivity 
to intervalley scattering are specific to all graphitic films 
independently of their morphology due to lower (trigo- 
nal) symmetry of the wavevector K in the corner of the 
hexagonal Brillouin zone of a honeycomb lattice crys- 
tal. Such a feature of the band structure does not af- 
fect the Josephson proximity effect in superconductor- 
graphene devices. This is because the propagation of 
a spin-singlet Cooper pair in graphene is related to the 
valley-symmetric Cooperon Cf , whereas Fermi statistics 
forbid the appearance of the other long-living mode C° 
in the superconducting order parameter. 
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